Higgs mechanism near the 5d bulk phase transition 



(N 



> 



(N 



Nikos Irges 1 , Francesco Knechtli 2 and Kyoko Yoneyama 5 

1. Department of Physics 
National Technical University of Athens 



O | Zografou Campus, GR-15780 Athens, Greece 

• 2. Department of Physics, Bergische Universitdt Wuppertal 

Q ■ Gaussstr. 20, D-^2119 Wuppertal, Germany 

(N 

e-mail: irges@mail.ntua.gr, knechtli@physik.uni-wuppertal.de, 
yoneyama@physik.uni- wuppertal. de 

43 



Abstract 



l/"") ! We present a non-perturbative model of Gauge-Higgs Unification. We consider a five- 

dimensional pure SU{2) gauge theory with orbifold boundary conditions along the fifth 



dimension, such that the symmetry is reduced to U(l) at the fixed points of the orbifold 
action. The spectrum on the four- dimensional boundary hyperplanes includes, apart from 
the U(l) gauge boson, also a complex scalar, interpreted as a simplified version of the 
Standard Model Higgs field. The gauge theory is defined on a Euclidean lattice which is 
anisotropic in the extra dimension. Using the boundary Wilson Loop and the observable 
that represents the scalar and in the context of an expansion in fluctuations around a 
Mean-Field background, we show that a) near the bulk phase transition the model tends 
to reduce dimensionally to a four-dimensional gauge-scalar theory, b) the boundary U(l) 
gauge symmetry breaks spontaneously due to the broken translational invariance along 
the fifth dimension, c) it is possible to construct renormalized trajectories on the phase 
diagram along which the Higgs mass is constant as the lattice spacing is varied, d) by 
taking a continuum limit in the regime where the anisotropy parameter is small, it is 
possible to predict the existence of a Z' state with a mass around 1 TeV. 
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1 Introduction 



In the Standard Model (SM) the Higgs field H is introduced as a fundamental scalar 
and inserted in the classical Lagrangean via the most general potential of engineering 
dimension four, consistent with the field's quantum numbers: 

V H = -fi 2 H^H + X(H^H) 2 . (1) 

The relative sign between the two terms is not fixed by any symmetry and is an external 
assumption. The Higgs mechanism then automatically proceeds by the field developing a 
vacuum expectation value (vev) < H >= v/y/2 which, upon minimization of the potential, 
turns out to be non-zero and satisfying v = /// y/2 at the classical level. At the same order, 
the Higgs mass is m H = v\/2\ and the neutral gauge boson mass m z = l/2vg with g 
a coupling derived from the gauge couplings of the group factors that contribute to the 
mass. Thus, in its minimal version, the Higgs mechanism is described by a gauge coupling 
g, a dimensionless quartic coupling A, a dimensionful mass parameter \x and the vev v. 
External input is also the potential itself, in the sense described above. Fluctuations 
around this non-trivial vacuum define the SM in its spontaneously broken phase, a state 
of matters that seems to be consistent even with the most recent LHC data. The origin of 
the ingredients that conspire to make the mechanism work is however unknown. Perhaps 
the most convincing clue that this cannot be the end of our formulation of the low energy 
end of high energy elementary particle physics is the quantum response of the fluctuations. 
There is a quadratic dependence of the scalar mass on the cut-off, believed to render the 
theory unnatural for a light Higgs particle. Historically the dominant solution to this 
puzzle has been supersymmetry. Here we propose an alternative scenario where the 
mechanism develops dynamically, described by three (in infinite lattice volume actually 
only by two) dimensionless parameters, without introducing an explicit potential or a vev. 
Moreover, we demonstrate that in our proposed scheme the mass of the Higgs particle is 
insensitive to the cut-off along renormalized trajectories, without supersymmetry. 

The general context is that of "Gauge-Higgs Unification" (GHU) [1, 2], where the 
Higgs field originates from the components of a higher dimensional gauge field along the 
extra dimension(s). Since we would like to have a control of the theory at the quantum 
level, we will exclude from our discussion warped and curved space-times. Instead, our 
starting point is a five- dimensional SU(2) gauge theory compactified on the S 1 /Z 2 orb- 
ifold, possibly the simplest prototype GHU model where the dynamical Higgs mechanism 
can be studied. The boundary conditions introduce four-dimensional boundaries at the 
fixed points of the orbifold, where the gauge symmetry is reduced to U(l). Subsequently, 
the boundary symmetry can, in principle, break spontaneously, generating a massive Z 
boson and a massive Higgs. Originally, this model was studied in the perturbative regime 
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where the first peculiar deviation from what we see in the SM was recognized to be the 
fact that fermions are necessary to trigger spontaneous symmetry breaking (SSB), tying 
the existence of SSB to the fermionic content [2, 3]. Another obscure fact that adds to 
the above is that even if one accepts this as a necessary feature of GHU models, it seems 
to be hard to have a Higgs heavier than the neutral Z gauge boson, without consider- 
ing values for the parameters involved somewhat unnatural. Also, in the perturbative 
approach, even though the scalar potential is a dynamical quantum ( Coleman- Weinberg- 
Hosotani (CWH)) effect, a vev must be inserted by hand, like in the SM. Finally, the 
non-renormalizability of the gauge theory often brings doubts about certain attractive 
conclusions drawn from perturbative loop calculations, such as the finiteness of the Higgs 
mass [4, 5]. From our point of view, it is this non-renormalizability along with the pertur- 
bative triviality of the five-dimensional gauge coupling (perhaps not an unrelated property 
to the former) that may be the root of these problems. 

For the above reasons, we have started a non-perturbative investigation of these simple 
orbifold gauge theories [6]. 1 At an early stage, an exploratory lattice Monte Carlo (MC) 
study of the £77(2) orbifold theory was performed that revealed that SSB is present 
already in the pure gauge system, signaled by a massive Z boson [12, 13]. At the same 
time however also the practical difficulty of a systematic MC study and the necessity for 
a non-perturbative analytic approach became obvious. Recently such a formalism was 
developed [14, 15], consisting of an expansion of the path integral in fluctuations around 
a Mean-Field (MF) background. Furthermore, introducing an anisotropy along the fifth 
dimension proved to be fruitful. There is increasing confidence by now that this MF 
expansion is a faithful representation of the non-perturbative system in five dimensions. 
If the five- dimensional SU(2) gauge theory has periodic boundary conditions along the 
extra dimension, dimensional reduction to four dimensions (if it occurs) drives the system 
to a Georgi-Glashow model (i.e. a Ad £77(2) gauge theory coupled to an adjoint scalar) 
from the lower dimensional point of view. Our Mean-Field analysis suggests that with 
these boundary conditions SSB does not occur. If instead we consider orbifold boundary 
conditions, due to the breaking of the £77(2) gauge symmetry to U(l) at the boundaries, 
dimensional reduction takes us to a four- dimensional Abelian-Higgs system. Our Mean- 
Field expansion method indicates that SSB is realized in this case [15]. In fact, the 
static potential along the boundaries is of a Ad Yukawa type and the (smallest) Yukawa 
mass corresponds to the mass of the Z gauge boson, confirming the earlier Monte Carlo 
simulations of the orbifold model. 

In order to understand the origin of SSB, we have to take a closer look at the struc- 
ture of our construction. On the left of Fig. 1 we show a schematic picture of the orbifold 

1 Recent Monte Carlo investigations of the periodic theory include [7, 8, 9, 10, 11]. 
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Figure 1: The lattice orbifold. The Mean-Field background depends on the extra- 
dimensional coordinate in a way that induces the spontaneous breaking of the gauge 
symmetry. Lighter shades of a given color indicate different MF background values 
v {n 5 + fe/2), k = 0, 1, 2, 3, • • • of the links. 

lattice we have in mind. Its precise definition has been given elsewhere [6]. The four- 
dimensional boundaries at the two fixed points of the orbifold action are covered by U(l) 
links (w-links). All other links are £77(2) links (s-links), except from those lying along 
the fifth dimension and whose one end touches the boundary (/t-links). These have one 
"end" transforming as U(l) and the other as SU(2). This is the proper orbifold lat- 
tice, invariant under the gauge and orbifold actions. In fact, all correlators representing 
physical observables like the Wilson Loop and the Higgs are also gauge invariant and in- 
variant under the orbifold action [6]. As a result, the breaking of translational invariance 
along the extra dimension is spontaneous. In a 1-loop CWH computation this breaking is 
encoded merely in a modification of the pre-factor of the potential, leading to the conclu- 
sion of no SSB, as on the torus. In the present approach on the other hand already the 
Mean-Field background contains the necessary information. The Mean-Field equations 
for the background produce a non-trivial profile for it along the extra dimension [16, 15], 
schematically represented on Fig. 1, where different shades of a given color represent dif- 
ferent MF background values of the links. On the right of the figure we have explicitly 
indicated a few of those values. The non-trivial profile is of course a direct consequence 
of the spontaneously broken translational invariance, allowing us to call sometimes the 



4 



MF background (together with the vibrations of the lattice above it) as a "phonon" in a 
condensed matter language. For similar reasons, the spontaneously broken translational 
invariance is felt by the system without the need for any background, in a full Monte 
Carlo study [12, 13]. 

Finally in regards to the connection of our MF approach to perturbation theory, we 
mention that each order of this expansion is believed to be equivalent to the summation 
of an infinite number of perturbative Feynman diagrams [17]. Practically, by trivializing 
the background and taking the lattice coupling to infinity, we can reach the perturbative 
regime and reproduce the corresponding results. 



2 The Higgs mechanism as a phonon trigger effect 

Our anisotropic lattice has L spatial points, T time-like points (with lattice spacing 04) 
and A^5 + 1 points along the fifth dimension (with lattice spacing 05). L will be removed 
by increasing it enough so that physics does not depend on it. Time will be used to 
extract the masses of the ground and first excited states. In short, the five-dimensional 
lattice coupling j3 = Aa^/gl (g 5 is the dimensionful 5d gauge coupling), the anisotropy 
parameter 7 = a 4 /a 5 (at the classical level) and N 5 are the three dimensionless parameters 
that parametrize our model. The action is the five-dimensional Wilson plaquette action, 
anisotropic in the fifth dimension. 

The quantities that will be used in our analysis are the Wilson Loop on the boundary at 
the origin and the Higgs observable. In [15] we computed them to first non-trivial order in 
the fluctuations around the MF background and we will not repeat the calculations here. 
We will instead compute these observables numerically and extract the static potential 
from the former and the Higgs mass from the latter. The specific quantities that we are 
interested in are 

F 1 = m H R, (2) 
the Higgs mass in units of the inverse size of the extra dimension 1/R, the ratio 

Phz = (3) 

m z 

of the Higgs to Z mass, with the Z mass extracted from a fit of the static potential to 
the form {b is a constant and r is the spatial length of the Wilson Loop) 

e -m z r 

V(r) = -b h const. (4) 

r 

and 

Phz> = , (5) 

m z > 
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the ratio of the Higgs mass to the mass of the first excited vector boson state. This 
is usually called a Z 7 , and it can also be derived straightforwardly from a fit of the 
static potential. Clearly, such a fit makes sense if the spectrum can be interpreted as 
an effective four-dimensional theory which by itself is not a precise enough definition 
of a satisfactory dimensional reduction. Our more constrained criteria for dimensional 
reduction are therefore that 

• the fit to Eq. (4) is possible with mz 7^ 0. This ensures that there is SSB, signaled 
by the presence of the massive U(l) gauge boson. Otherwise the gauge boson is 
massless and only a Coulomb fit is possible. 

• the quantities Mh = a^mu and Mz = a^mz are < 1. This ensures that we are 
in a regime of the phase diagram where the lattice spacing does not dominate the 
observables. 

• we have muR < 1 and pnz > 1- These two conditions ensure that the Higgs and the 
Z mass are lighter than the Kaluza-Klein scale 1/R on one hand and that the Higgs 
is heavier than the Z on the other, a desirable situation from the phenomenological 
point of view. In fact, we will target the value 

Phz = 1-38, (6) 

which is (approximately) the currently favored value of the analogous quantity in 
the SM, based on recent LHC data [18, 19]. 

In summary, we have the observables iq, p HZ and pnz r , all three depending on the three 
dimensionless parameters (3, 7 and N 5 . Our method then is to fix iq to a given value 
and keep puz fixed to the value Eq. (6). This leaves puz 1 be a function of one parameter 
which we choose to be N 5 and by doing so we obtain a value for the mass of the Z' for each 
N 5 . We call such a trajectory on the phase diagram that also fulfills our three conditions 
for dimensional reduction, a Line of Constant Physics (LCP). For illustrative purposes we 
will also give dimensionful values for the masses by inserting in the ratios the SM value 
for m z . 

Eventually we would like to understand the physical meaning of the N 5 — > 00 limit and 
for that we have to describe the structure of the phase diagram. In the MF expansion the 
phase diagram can be plotted already at the level of the background. Even though one 
could consider corrections due to fluctuations, we will stay at the lowest order, because 
the corrections can be seen to be small. The background value of the gauge link variables 
on the anisotropic lattice are denoted as vo(n$) along the p = 0,1,2,3 directions and 
as v 05 (n 5 + 1/2) along the fifth dimension, with n 5 denoting the corresponding integer 
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coordinate. The phases of the system can be defined from the phonon profile as follows 
(the statements hold Vn 5 ): 



• Confined phase: v (n 5 ),v 05 (n 5 + 1/2) = 0. 

• Layered phase: Vo(n 5 ) 7^ 0, ^05(^5 + 1/2) = 0. 

• Deconfined phase: ^0(^5)^05(^5 + 1/2) 7^ 0. 



According to the MF method, the boundary between the Deconfined phase and the Con- 
fined and Layered phases has a different order depending on (3 and 7. Tuning (3 so that 
one is in the Deconfined phase and always near the phase boundary, one finds that for 7 
larger than a value that is slightly less than 1, the phase transition is of first order, while 
below that value it turns into second order. We emphasize that the phase transitions 
defined in this way are always bulk phase transitions, a fact that has been extensively 
verified on the periodic lattice [10]. Even though it is much harder to verify via Monte 
Carlo simulations the change of the order of the phase transition at small 7 [9], we will 
assume here that the order of the phase transitions that the MF predicts is always correct. 

When the bulk phase transition is of first order, the four-dimensional lattice spacing 
04 remains finite no matter how closely the phase boundary is approached. Then, since 
(/ is the physical four- dimensional volume) 

L = —, N 5 = — 7 , (7) 

04 04 

in the L = N 5 = 00 limit the physical size of the system goes to infinity at a finite lattice 
spacing. When the phase transition is second order, one expects instead that the lattice 
spacing goes to zero at a finite physical volume. There is an important physical difference 
between the two cases. In the first case the low energy theory is an effective theory that 
must be defined with a finite cut-off. The existence of the LCP nevertheless guarantees 
that a sensible (in the quantum sense) and predictive effective theory exists. In the second 
case SSB would persist in the continuum limit and there is no need for a cut-off in the 
effective action. This would render the theory non-perturbatively renormalizable. 

With or without continuum limit, we will show that in the vicinity of the phase bound- 
ary we have a gauge-scalar system without a hierarchy problem and with a dynamical 
Higgs mechanism that may be described as a phonon trigger effect: as in crystals where 
the formation of Cooper pairs takes place and whose interaction with the phonon leads 
to superconductivity, here the Polyakov Loop with quantum numbers appropriate for the 
Higgs boson (of charge 2 under the U(l), see [13]) interacts with the MF background 
triggering the spontaneous breaking of the gauge symmetry. Like on the periodic lattice, 
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Fx 


N 5 


7* 


P* 


0.61 


12 


0.5460(33) 


1.343501425 




14 


0.5320(10) 


1.34442190 




16 


0.5228(7) 


1.34664820 




20 


0.5028(18) 


1.35582290 




24 


0.4844(32) 


1.36940695 


0.20 


6 


0.5113(15) 


1.351160631 



Table 1: Bare parameters of the LCP defined by pnz = mu/mz = 1.38 and Fi = uihR = 
0.61, together with one point for a LCP with puz = 1-38 and F\ = 0.20. The lattice gauge 
couplings j3* correspond to the central values 7* and are computed for future reference. 

our motivation to believe that the lattice spacing decreases as the phase transition is ap- 
proached is that both Mu and Mz decrease. In fact, the only place where one can have 
Mh, Mz < 1 is reasonably near the phase transition. In the small 7 regime one can have 
values pretty much as small as desired. In the large 7 regime on the other hand there 
seems to be a barrier below which the masses can not be decreased anymore. Therefore 
it is not a surprise that if it is at all possible to construct an LCP as defined above, it will 
surely be a line near the phase boundary. 

3 Lines of Constant Physics and the Z' 

The first LCP we construct is one where 

Fi = m H R = 0.61 , p HZ = 1-38 (8) 

are kept fixed. Along this LCP, we compute p HZ > for N 5 = 12, 14, 16, 20, 24. On Fig. 2 we 
plot the corresponding points interpolated by a black line on the phase diagram, which 
are listed in Table 1. As discussed, it is a line near the phase boundary and in particular 
in the small 7 regime where the phase transition is of second order according to the MF. 
Thus we can attempt to take the continuum limit. Note that this is also the regime 
where from both MF [20] and MC [10] studies of the periodic lattice it has been seen that 
gauge fields are localized on four-dimensional hyperplanes. Evidently this is an effect 
independent of boundary conditions. 

For each value of jV 5 , we compute the Z and Z' masses for various values of the 
parameter 7. The third parameter (3 = j3(j, N 5 ) is set by requiring that i*\ has the 
desired value 0.61. The gauge boson masses are extracted by identifying them as Yukawa 
masses describing the static potential V(r) along the boundary. From the static force 
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P y 

Figure 2: LCP (black line) defined in Eq. (8) near the (tricritical point of the) bulk phase 
transition. Red: Confined phase. Blue: Layered phase. White: Deconfined phase. The 
magenta point (star) is on a different LCP with F\ = 0.2, pnz = 1-38. 

F = dV(r)/dr we compute the quantity 

y '(r) = [y(r + a 4 ) - y(r - o 4 )]/(2a 4 ) , y(r) = \n(r 2 F(r)) . (9) 

A Yukawa mass is identified as a plateau in the quantity —y'{r) as a function of r. As it is 
shown on the left plot in Fig. 3 for N 5 = 24 and 7 = 0.485, typically there are two plateaus, 
the smaller (red points) we take to define Mz and the larger (blue points) Mz>- The ranges 
of r values defining the plateaus are taken around the minima of the derivative of —y'(r). 
The masses are the average and the errors are the standard deviation of the plateau points. 
The value of Mz is improved by iteratively solving Mz = — [o,^y'{r) — Mz/(Mzr/a^ + 1)]. 
The value of L should be large enough to clearly identify the plateaus and we set L = 400 
for all N$ values. 

In the next step, at fixed iVs, we compute puz as a function of 7 from the known values 
of Mz and Mjj. The data can be very well fitted by a straight line, as is demonstrated 
for A^5 = 24 by the upper red points and the red line on the right plot of Fig. 3. From 
the linear fit we compute the value 7*(A 5 ) which gives the desired value puz = 1-38. The 
error on 7* takes the correlation of the fit parameters into account. For the data shown in 
Fig. 3 we get 7* (24) = 0.4844(32). A summary of the LCP parameters for all N§ values 
is given in Table 1. 
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Figure 3: Left plot: plateaus of the quantity defined in Eq. (9) corresponding to the Z 
(red points) and Z' (blue points) masses. Right plot: the puz data (upper red circles) are 
lineraly interpolated (red line) to the value of 7 corresponding to puz = 1-38 (marked by 
the dashed horizontal line). The lower blue circles show the data for pnz 1 with a linear 
fit (blue line). 

Finally we compute, for each value of N 5 , the ratio puz 1 when 7 is set to the value 
7*. This is done by fitting linearly in 7 the data of pnz 1 as shown by the lower blue 
points and blue line on the right plot of Fig. 3. We take the fit result evaluated at 7*(A^" 5 ) 
and augment the fit error by adding the slope of the fit multiplied by the uncertainty 
of 7*(iV 5 ). For the data shown in Fig. 3 we get Phz>{i*) = 0.2313(1)(34), where the 
dominant error comes from the uncertainty of 7*(iV 5 ). On Fig. 4 we plot our data of pnz' 
against a^mu for N 5 = 12, 14, 16, 20, 24. Notice that since F\ = mnR = (a>4jnH)N$/ (t*7t), 
atfriH is proportional to j*/N 5 along the LCP and it is a 4 m# which measures the physical 
distance to the continuum limit (it is the inverse correlation length). In principle we could 
fit the points with a quadratic curve because from the Symanzik analysis of cut-off effects 
we expect the dominant contributions to come from the dimension 5 boundary operator 
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(10) 



multiplied by one power of the lattice spacing and from the dimension 7 bulk operator 

h h E tr {^mn (D 2 m + D%) F MN ) (11) 

^ 5 M,N 

multiplied by two powers of the lattice spacing [21]. In fact, because we are very close to 
the phase transition, we are in a regime where the effect of the dimension five boundary 
operator dominates, thus the linear fit on Fig. 4. The extrapolation to a^rrtu — > leaves 
a non-zero intercept with the vertical axis which corresponds to 



p HZ > = 0.1272 



(12) 
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p HZ =1 .38, F 1= 0.61 , on the boundary 




0.05 



0.02 0.04 0.06 0.08 0.1 

a 4 m H 

Figure 4: Extrapolation of LCP in Eq. (8) to a^rriH — > 0. 

For mz = 91.19 GeV this implies a Z' of mass mz> = 989 GeV in the continuum limit. 
The x 2 P er degree of freedom of the fit is excellent, 0.025/3. 

Before we state our conclusions, it is worth making a comparison between the SM Higgs 
mechanism and our scheme. As these are two different theories, there is no well defined 
way to do this so we can only be qualitative. The dimensionless 4d gauge coupling g can 
be thought to be analogous to the dimensionless 5d gauge coupling (3. The Higgs field H 
is introduced in the SM by hand while in the lattice GHU model its presence is induced 
by the orbifold boundary conditions, which reduce the adjoint scalar to the complex 

1 2 m 

scalar made from A 5 ' . The quartic coupling A of the SM Higgs sector is a dimensionless 
parameter like the anisotropy parameter 7. The choice of the negative sign in front of the 
fi 2 term in the Higgs potential triggers SSB in the SM. In the MF expansion we see SSB 
once the background is defined as the point around which the path integral is expanded, 
a signature of the broken translational invariance in the extra dimension. The presence of 
the phonon seems to be the crucial fact that triggers the subsequent spontaneous breaking 
of the gauge symmetry, like in superconductors. In the SM it is likely that the parameter 
Phz takes the approximate value 1.38, an experimental fact. On our lattice orbifold there 
is a family of pnz = 1-38 LCP's, each labelled by a different value of tuhR and we have 
chosen to plot one of these, the one that corresponds to F\ = muR = 0.61. We have 
checked that there is in principle no obstruction in constructing an LCP with p#^=1.38 
and a much smaller F\ (we have generated a point for F\ = 0.20, see Table 1 and Fig. 2 
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Figure 5: Same as in Fig. 3 but for Fi = 0.20. 



and Fig. 5), but numerically this is slightly more demanding since the masses in lattice 
units become very small as a function of N$. 

We checked also the region where 7 is one or larger and found that an LCP with 
Puz — 1-38 is possible only in the small gamma regime. In the regime 7 > 1 dimensional 
reduction can occur only through compactification and therefore we choose F\ between 
0.10 and 0.20. For 7 around one we get p^z values which are smaller than but close to 
one, in agreement with the results from Monte Carlo simulations at 7 = 1 [12, 13, 21, 22]. 
The pnz' values are ~ 0.090. For 7 = 4 we see only one plateau for the quantity defined 
in Eq. (9), which we interpret as the Z mass. The puz values are again consistent with 
one. The fact that we do not see a second plateau for the Z ! mass may be reasonable 
since we are in the compact phase and it could be that the Z ! state is too heavy. 



4 Conclusions 

We presented a numerical analysis of observables computed in an analytical expansion in 
fluctuations around the Mean-Field background on the anisotropic lattice orbifold, defined 
and developed in previous publications. We showed that non-perturbatively spontaneous 
symmetry breaking is a property of the pure gauge system and we were able to draw 
lines on the phase diagram along which the ratio of the Higgs over the Z boson mass 
puz remains finite. Furthermore, by taking the continuum limit near the bulk phase 
transition and at values of the anisotropy parameter around 7 ~ 0.5, we demonstrated 
that the first excited state in the vector boson sector is light, with a mass in the TeV 
regime. The Standard Model value of puz — 1-38 for which this Line of Constant Physics 
was drawn, could be reached only in this, small 7 regime. Even though we used a toy 
SU (2) model, we believe that these generic properties will persist in more realistic cases 
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where the bulk gauge symmetry may be for example 577(3) x SU(3) or 577(5). Similarly, 
we expect that the presence of fermions will not alter the observed qualitative behavior. 
We plan to perform these generalizations in the future. The scheme of the Mean-Field 
expansion seems to be a good (semi) analytical description of the non-perturbative system 
in five dimensions thus it should be considered at least as a valuable complementary tool 
in the study of Gauge-Higgs Unification. For among others, it provides us with a guide 
to a Monte Carlo study, pointing to the regime on the phase diagram where one should 
perhaps focus. Especially if the absence of low energy supersymmetry is experimentally 
confirmed, our results here may have given us a hint for an alternative solution to the 
Higgs hierarchy problem and for a possible dynamical, non-perturbative origin of the 
Higgs mechanism. 
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